Slow relaxation and aging phenomena at nano-scale in granular materials 
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Complex systems having metastable elements often demonstrate nearly log-time relaxations and 
a kind of aging: repeated stimuli weaken the system's relaxational response. Granular matter is 
known to exhibit a wealth of such behaviors, for which the role of thermal fluctuations is usually 
ignored. However, we demonstrate that the latter can pronouncedly affect contacting mesoscopic- 
scale asperities and be macroscopically observed via appropriate acoustic effects. We also propose a 
mechanism comprising slow relaxations and aging as intrinsic properties of a wide class of systems 
with metastable states. 

PACS numbers: 83.80.Fg, 83.10.Tv, 45.70.-n, 43.25.+y 



Introduction - Slow relaxation phenomena in granu- 
lar systems attract much attention due to both high im- 
portance of granular materials for geophysics and many 
industries and due to general physical interest to be- 
havior of complex, glassy-type systems for which gran- 
ular materials act as macroscopic analogs of ensembles 
of individual atoms and molecules [l|. In particular, 
slow compaction of granular materials due to mechan- 
ical perturbations has been widely studied and it 
has been shown that, in such a process, the response 
of the material is related to -sporadic- grain rearrange- 
ments. Compared to slow macroscopic {i.e. grain-scale) 
rearrangements in granular materials, processes on the 
mesoscopic level {i.e. below micrometer, down to nano- 
scale) driven by natural thermal fluctuations are much 
less studied. The latter are reasonably believed irrel- 
evant to grain rearrangements during compaction and 
jamming- unjamming transitions [l|, but using ap- 
propriate acoustic techniques, spontaneous thermally- 
activated mesoscopic-scale processes can also be quite 
macroscopically observed. In particular, observations of 
slow relaxation of the elastic- modulus value in laboratory 
samples with cemented granular structure [|| as well as 
similar effects in field measurements in sandy soil on a 
scale of ~ 10 m 9[ are known. High-intensity acoustic 
"conditioning" j8| or a mechanical impact [t| perturbed 
weakest bonds and produced perturbations in the mate- 
rial elastic modulus of order 10 -6 — 10~ 3 that were rather 
problematic to monitor. To overcome this difficulty, it is 
desirable to reveal a parameter that is dominated not by 
the stable material skeleton, but precisely by the per- 
turbed weak-bond-network. To this end, we make use of 
an acoustic-signal component produced by the own ma- 
terial nonlinearity, which is strongly dominated by the 
contributions of the weakest-contact fraction. We show 
that aging of granular materials induced by mechanical 
perturbations can be explained through a physical mech- 
anism, based on the bistable character of contacts. 



Methods - Feasibility of this idea was demonstrated 
in [l(| by using the so-called nonlinear cross-modulation 
technique to monitor structural perturbations in granular 
material bulk induced by weak mechanical shocks. An- 
other, practically simpler, nonlinear-demodulation tech- 
nique was successfully applied for studying fine structural 
changes - avalanche precursors - in slowly tilted granular 
packings ll|,|l2|. Here, the sounding technique 11 Jl2 | 



combined with pulse-type perturbations like in [lCj is 
used to study slow relaxation of the weak-bond network 
in gran ular material with particular attention to revealed 
in [12| aging of material undergoing repeated perturba- 
tions. The experimental system is a 1 mm or 2 mm glass 
bead packing in a container to which a small electro- 
magnetic shaker is attached. It produces rather weak 
(comparing to typical conditions of tap- induced com- 
paction [lj-ljl) perturbing pulses. The strain amplitude of 
the perturbing pulses is varied in different measurements 
from about 10 -7 to 1CP 6 and their duration is of 20 ms. 
The primary AM-modulated wave is at strain amplitude 
£a ~ 10 -8 — 10 -7 . Unlike observations of the primary- 
wave (used in [1, H[) dominated by the medium skeleton, 
we use the demodulated component that is strongly dom- 
inated by the weakest contacts in the material, because 
nonlinearity of the weakest, least- loaded Hertzian con- 
tacts strongly exceeds nonlinearity of stronger contacts 
forming the stable skeleton of the material [lg, [l7 1 . 



FigurcQ] shows various examples of slow relaxation 
in granular materials, observed using their nonlinear- 
acoustic response. As a function of time, the latter 
demonstrates power-law rates close to l/t 1+n with \n\ <C 
1, i.e., close to log-time behavior corresponding to n = 0. 
Plots (a) and (c) demonstrate peculiar weakening of the 
material reaction to series of identical taps, i.e., a kind of 
"aging" . Plots (a) and (b) also show much higher sensi- 
tivity of the nonlinearity-produced signal to the state of 
broken and then gradually restored weakest bonds com- 
pared with the fundamental component variability. As 
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FIG. 1. Various examples of nonlinear-acoustic monitoring 
of slow relaxations: (a) - after perturbing shocks using cross- 
modulation technique [l(| ; (b) - due to internal micro-slips in 
gradually tilted granular packings observed using demodula- 
tion technique [H], Q3i ( c ) ~ similarly observed "aging" of the 
system response to a series of identical fairly weak perturb- 
ing shocks, and (d) - relaxation after one stronger shock and 
after different-amplitude shocks that perturb barriers with 
essentially different energies. Much weaker perturbation of 
fundamental (linear) component is shown in plots (a) and (b) 
for comparison. 
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FIG. 2. (a) Adhesion/compression force (solid line) for a 
small asperity (r = 30 nm with elastic parameters and surface 
energy typical of silica and glasses) and equilibrating elastic 
force of a cantilever (dashed lines) which jointly are able to 
produce bistable equilibria, (b) Elucidation that Hertzian-like 
elastic tensile force for sufficiently large contacts subjected 
to short-range adhesion forces from an opposite solid sur- 
face can also form bistable equilibria, (c) Corresponding two- 
minimum potential wells, for which the energy-barrier depths 
Eb, c exhibit opposite trends as a function of separation A nor- 
malized to atomic size d a - (d) Illustration that those energy 
barriers should be overcome to break "closed" contacts and 
close "open" ones. Dashed rectangles show contacts close to 
bistability-zone boundaries that can be sensitive to thermal 
fluctuations. 



argued in [ll|, such rapid and strong drops in the de- 
modulated signal are related to breaking of the weak- 
est contacts. Their nonlincarity is significantly higher 
than that of stronger contacts composing the material 
skeleton and dominating its linear elastic moduli. Thus 
the post-shock signal increase reflects restoration of bro- 
ken contacts. Note that even if the probing signal is 
switched off just after the shock and switched on after a 
pause, the nonlinearity restores spontaneously, showing 
that the probing wave influence does not dominate the 
effect, although high-intensity acoustic strains (say 10 -5 ) 
may even be able to perturb the weak bonds as in [1, [l]| . 

Mechanism - Under room temperature T ~ 300 K the 
characteristic thermal energy fcgT (ks is the Boltzmann 
constant) unambiguously indicates that thermal fluctu- 
ations cannot affect the state of visible (even weakly 
loaded) macroscopic contacts usually considered [H, 0] 
in granular matter modeling. Thus only mcsoscopic sur- 
face asperities (from tens to hundreds of nanometers) can 
be considered as candidates of bistable structural ele- 
ments potentially sensitive to thermal fluctuations. To 
understand the origin of their bistability, the analogy 
with the bistable behavior of tips in atomic-force mi- 
croscopy (AFM) is very useful. For a tip already com- 
pressed by the contacting solid, the elastic Hertzian force 
is repulsive, whereas the tip yet approaching a solid sur- 
face experiences the influence of short-range attraction 
forces (often modeled using the Lcnard- Jones potential). 
This attraction force for a AFM tip approaching another 



solid is equilibrated by the elasticity of the cantilever 
(dashed lines in Fig. 2a). If the cantilever is soft enough, 
a bistablity zone can appear [l4[ as illustrated in Fig. 2a. 
In this zone, for a given initial position A\ < A < A2 of 
the unstressed cantilever, it can equilibrate the attraction 
force at two positions of the tip, "closed" and "open" , (in 
the latter, the attraction force is almost absent). If the 
cantilever is moved forth and back, peculiar hysteretic 
jumps between the two positions occur (see arrows in 
Fig. 2a). Inside the bistability zone eventual perturba- 
tions (e.g., sufficiently strong thermal fluctuations) can 
cause transitions between the two equilibrium states. 

At first glance, for mesoscopic asperities at grain sur- 
faces, there is no "soft cantilever" to create bistable equi- 
libria by analogy with AFM. However, one can recollect 
general arguments jl5| leading to the fact that for com- 
pressed contacts (not necessarily ideally spherical), the 
stored elastic energy scales as h 5 ^ 2 (where h is the dis- 
placement of the contact apex). This means that for a 
compressed contact, the elastic force F comp scales by the 
famous Hertzian law F comp ex ft 3 / 2 . But the same argu- 
ments applied to a contact apex displaced by the value 
\h\ = \x — A\ due to a localized attractive (stretching) 
force also lead to appearance of a tensile elastic force 
Ftens oc \x — A\ 3 / 2 that equilibrates the attraction (see 
dash-dotted curves in Fig. 2b). 

Unlike AFM cantilever elasticity, this elastic force is 
nonlinear, i.e., initially it can be sufficiently soft to ere- 
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FIG. 3. AFM image of micro-asperities (a) and schematically 
shown action (b) of perturbing shocks and post-shock relax- 
ation in the diagram of energy barriers Eb, c - The lengths of 
vertical arrows depict different transition rates. The boxes 
depict active energy ranges where the contact population sig- 
nificantly changes due to thermally-assisted transitions. 



ate the second (distant) equilibrium position for the con- 



tact tip. For AFM tips with typical radius 



< 



10 nm, 



however, this "nonlinear spring" hidden inside the tip 
is insufficiently soft relatively to adhesion, so that only 
in-sequence connected soft cantilever can create the sec- 
ond potential minimum. But for a larger contact radius 
r > 30 — 50 nm, due to a different dependence on r for the 
attraction and the "hidden spring", the latter becomes 
sufficiently soft relative to the adhesion force (compare 
the curves for r = 30nm and r = 200nm and other iden- 
tical parameters in Fig. 2(b)). Thus, bistable equilibria 
with two-minima potential wells can appear in a finite 
range of separations between the asperity tip and the 
opposite surface without a real soft cantilever. Fig. 2(c) 
schematically shows how the resulting two-minima po- 
tential evolves with the initial separation A. These rep- 
resentations suggest physically clear interpretation to the 
well-known 18( transition from the so-called DMT-model 
of very small contacts (that do not exhibit adhesion hys- 
teresis) to JKR-modcl for larger contacts exhibiting ad- 
hesion hysteresis. The latter can be viewed as a special 
case of mechanical hysteresis (much like for AFM tips, 
but without the necessity of an artificial soft cantilever) . 

Assuming elastic and surface energy constants typical 
of glass-like materials, one concludes that for contacts 
with r ~ 10 2 nm and room temperature T ~ 300K, there 
are narrow regions near the boundaries of the bistability 
zone, where one of the potential wells Ej, or E c (Fig. 2b) 
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FIG. 4. Evolution of the number of closed mesoscopic contacts 
(left) and populations Nb(E) of broken ones (right). Labels 1 
are for the initial state of just prepared packing with Nb(E) = 
1 relaxed during 30 s ((a) and (b)) and during 1500 s ((c) and 
(d)). Labels 2 are for the system state just after a fairly weak 
perturbing shock that breaks only a portion of earlier closed 
contacts. Labels 3 correspond to the system states 24s after 
the shock. Compare experimental Fig. Id with Fig. 4a, where 
the dashed curve shows continuous relaxation. 



firmed the presence of numerous asperities about 10 2 nm 
in radius and 20— 50nm in height (Fig. 3a) consistent with 
the values reported in . A single contact between two 
grains leads to 10 3 — 10 4 micro-asperities. Even if 1% 
of them actually gets in contact, one obtains ~ 10 2 of 
such mesoscopic contacts for a visible one. Following 
we conclude that contribution of such loose but numer- 
ous micro-contacts can dominate over the nonlinearity of 
much stronger (and thus less nonlinear) macro-contacts 
creating the material skeleton. This explains why non- 
linearity can drop drastically after fairly weak perturba- 
tions that still leave the material skeleton intact, but can 
suffice to break the mesoscopic contacts. 

Figure 3b shows the physical meaning of relaxational 
closing of open contacts and the destructive action of 
perturbing weak tensile shocks (which do not completely 
get the system out of the bistability region). We recall 



is much greater than kgT (say 10 — 10 times), whereas that even for large m icro-contacts (with r > 10 2 nm) 



the other one is ~ 10 1 fcsT (see Fig. 2d for r = lOOnm) 
Thus, near the left boundary of the bistability region 
the closed state is much stabler, whereas the open one 
is metastable. Near the right boundary, the situation 
is opposite. The metastable equilibrium energy is com- 
parable with that of thermal fluctuations, so that the 
latter are able to induce jumps to the opposite stabler 
state with characteristic waiting times tq expiE^c/ksT) 
according to the Arrhenius law, where the attempt time 
To for nanometer-scale tips of the asperities can reason- 
ably be taken tq ~ 10~ 12 s. 

Direct AFM inspection of the glass-bead surfaces con- 



the bistality zone is of the order of atomic size. Note 
that for characteristic attempt times tq ~ 10~ 12 s, only 
"active" contacts with barriers < can participate 

in thermally-induced transitions, for which the waiting 
times are below tens of hours. Then the narrow "ac- 
tive" parts of the energy curves near the bistability- 
region boundaries can be reasonably well approximated 
by straight segments (thick solid lines in Fig. 2d) . Con- 
sequently, in such narrow regions for almost arbitrary 
distributions of asperities' heights, the density of energy 
states for the active mesoscopic contacts can be approx- 
imated as constant. 
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Kinetic Monte Carlo - We simulated the above dis- 
cussed transitions between "open" and "closed" contacts 
(their total number in simulations is 3 • 10 4 ) using a ki- 
netic Monte-Carlo approach. It uses Arrhenius transi- 
tion probabilities characterized by the energy barriers 
E c for jumps into closed states and Eb to broken states. 
To match experimental time scale we restrict ourselves 
to < Eb, c 45/sbT. If initially all mesoscopic con- 
tacts are broken, the population of closed contacts is 
zero, N C (E) = 0, whereas the broken states are fully 
populated, Nb(E) = 1. Gradual closing (relaxation) of 
the broken contacts looks as the motion of the "clos- 
ing front" of Nb{E) from small energy barriers to larger 
ones towards the right boundary of the bistability zone 
(as shown in Fig. 4b and 4d, curves 1, for two moments 
30s and 1500s). The rate of the front motion exponen- 
tially decreases with time. We recall that the energy 
scale of the abscissa is proportional to distances from 
the bistability-zone boundaries. This means that tensile 
perturbations temporarily shift the micro-contacts (in- 
cluding those corresponding to the current position of the 
"relaxation front" ) towards the right boundary of the en- 
ergy diagram. There, the energy wells for the open states 
are deep, whereas those for closed states are shallow (see 
Figs. 2c and 3b), so that transitions of previously closed 
contacts back to open states are fostered. 

As a result, the pre-shock position of the closing front 
(Fig. 4b, d, curves 1) by the ending of the shock becomes 
shifted back to the left (Fig. 4b, d, curves 2), but then the 
closing front continues its movement to the right. For 
the same shock amplitude and duration, the resulting 
amount of the broken contacts that increase the popula- 
tion Nb(E) essentially depends on the position of the clos- 
ing front just before the shock. Notice the much smaller 
difference between curves 2 and 1 (Figs. 4b and 4d) for 
the initial relaxation time 30s than for 1500s. Sufficiently 
strong perturbations simply shift all bistable contacts to 
the right beyond the bistability zone (where only open 
states exist) and break all earlier closed contacts. Since 
displacements of contacts in granular materials for grains 
of radius R can be estimated as sR for mechanical strain 
e and since the width of the bistability zone is of the or- 
der of an atomic size, such "strong" shocks correspond to 
strains > 10~ 6 — 10 -5 for millimctric grains. For weaker 
shocks, the system response can be rather non-trivial and 
multi-variant depending on shock amplitude, duration, 
and previous history, e.g., the number of previous shocks 
and inter-shock intervals. For example, besides the differ- 
ence in perturbation strengths, Fig. 4 shows that for the 
same post-shock relaxation time, the closing front 3 gets 
already to the right from its initial position 1 in panel 4b 
in contrast to the opposite situation in panel 4d. 

Next, it should be recalled that besides the difference 
in the barrier energies Eb >c the widths of bistability re- 
gions can strongly differ for different-size contacts. Thus, 
the same perturbation can be "strong" for smaller con- 
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FIG. 5. Simulated example of "aging" of a well relaxed (dur- 
ing 3000s) system subjected to a series of weak perturbations. 
Plot (a) is the time-domain representation similar to the ex- 
perimental curve Fig. lc; (b) and (c) are the energy spectra 
of broken contacts. The corresponding moments in plots (a), 
(b), and (c) are numbered. The dashed lines demonstrate 
that initial perturbation N£(E) — N£(E) is much stronger 
than Ni!(E) — N^°(E) in the "aged" state, where the initial 
spectrum N^(E) almost coincides with the final N^ 1 (E) . 



tacts and "weak" for larger ones (as defined above), so 
that perturbation/relaxation regimes for such fractions 
of bistable elements arc quite different. Similar situa- 
tions can occur, e.g., for magnetic domains and other 
mesoscopic structural units that are not identical unlike 
individual atoms (molecules) of the same type. 

For the above-discussed features, the "aging" of relax- 
ation response to fairly weak perturbations is a natural 
consequence. Namely, a series of periodical moderate- 
strength shocks applied to previously well-relaxed ma- 
terial (Fig. 4c and d) gradually shifts the system state 
towards the one shown in Fig. 4a and b. Nevertheless, 
the regime shown in Fig. 4a and b (for which the re- 
laxed front 3 lies to right from the initial position 1) is 
not reached, because the system response saturates when 
the relaxation between the shocks becomes able to heal 
the perturbation ANb produced by every previous shock. 
This saturated value ANb caused by repeated weak per- 
turbing action of a fixed amplitude becomes significanly 
smaller than the reaction ANb of the well-relaxed system 
to the initial shocks. The transition to such "aged" re- 
action of the system is shown in Fig. 5 for the simplest 
situation of identical mesoscopic contacts for which the 
barriers differ only due to different initial separations A) . 
This model already fairly well reproduces the experimen- 
tally observed (see Fig. lc) gradual "aging" of the system 
response to repeated weak shocks. 

The discussed mechanism of micro-contact destruc- 
tion/restoration does not necessarily imply impact-like 
perturbations. It is directly relevant to the observations 
of avalanche precursors in slowly tilted granular pack- 
ings 
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[12| that can be viewed as internal "shocks", 
after which the relaxation is also due to the discussed 
mechanism. Furthermore, sufficiently gentle tilt-induced 
strain variations can also cause slow breaking of the meso- 
scopic contacts as shown in Fig. lb. If the critical angle 
is not reached and forth-and-back tilts are periodically 
repeated (with ~ lmin. period), tilt- induced variations 
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in the demodulatcd-signal amplitude strongly decrease 
(become "aged"). However, if the system is again tilted 
after ~ 1 hour rest (during which the broken mesoscopic 
contacts spontaneosly restore), the excursion of the sig- 
nal variations also significantly restore (see [13). Very 
high sensitivity of this technique is related to the fact 
that atomic-scale (and even sub-atomic) displacements 
control the observed nonlinear-acoustic response. 

In the above discussion of the nonlinear-acoustic re- 
sponse determined by cumulated contribution of large 
contact ensembles we focused on weak perturbations, 
after which the relaxation of the ensemble was consid- 
ered as an independent act for each mesoscopic contact. 
However, estimates show that if a significant portion of 
such contacts is broken, the surfaces of the macroscopic 
inter-granular contacts can experience atomic scale sepa- 
ration (comparable with or even greater than individual 
bistability-zone widths). This means that for stronger 
perturbations, the jumps of mesoscopic contacts between 
the bistable equilibria are not independent, and impor- 
tant role should be played by the secondary, really col- 
lective hysteretic/relaxational mechanisms that will be 
discussed elsewhere. Nevertheless, even for the small- 
perturbation regime, the revealed mechanism of indepen- 
dent nearly-logarithmic (or power-law) relaxations en- 
sures a rich variety of possible regimes. In particular, 
the aging effect has a rather general nature and should 
manifest itself for a wide variety of complex systems con- 
taining large ensembles of relaxing mesoscopic units. 
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